PINNING OF INTERFACES IN A RANDOM ELASTIC 
MEDIUM AND LOGARITHMIC LATTICE 
EMBEDDINGS IN PERCOLATION 
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Abstract. For a model of a driven interface in an elastic medium 
with random obstacles we prove existence of a stationary positive 
supersolution at non- vanishing driving force. This shows the emer- 
gence of a rate independent hysteresis through the interaction of 
the interface with the obstacles, despite a linear (force=velocity) 
microscopic kinetic relation. We also prove a percolation result, 
namely the possibility to embed the graph of an only logarithmi- 
cally growing function in a next-nearest neighbor site-percolation 
cluster at a non-trivial percolation threshold. 



I. Introduction and the main result 

In this article, we consider a model for the propagation of one- 
dimensional fronts immersed in an elastic medium subject to an ex- 
ternal driving force and randomly distributed obstacles. The goal is 
to understand the overall macroscopic behavior of such fronts and its 
dependence on the external forcing. Here we prove existence of sta- 
tionary solutions at positive driving force and thus the emergence of 
hysteresis. 

In order to precisely state our model, let (Q, B, P) be a probability 
space, u G fi. The random front at time t is given as the graph 
{x,u(x,t,oj)) of a function u : IRx(0,oo)xf2 — > R solving the semilinear 
fractional diffusion problem 

(1) u t (x, t, uj) = — (— A) 1 / 2 -u(x, t, u) — f(x, u(x, t,U)),0j) + F 
u(x, 0, uj) = 0. 

The function f(x,y,uu) > is assumed to be locally smooth in x and 
y for any u> and of the form of localized obstacles of identical shape 
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Figure 1 . Pulling sandpaper out of a glass of water. 



and random positions with uniform density, i.e., the obstacle centers 



are given by a 2-dimensional Poisson process. See Assumption |1.1| for 
a precise statement. The constant term F is an external loading, the 
fractional Laplacian models the interaction of the front with the elastic 
medium in which it is immersed. 

Evolution problems of this kind arise in a large number of physical 
systems. A particularly simple example is that of "pulling sandpaper 
out of a glass of water". As illustrated in Figure [TJ we model the 
evolution of the wetting line of the water surface on a rough plate as 
the plate gets pulled out of the water. Equation (JT]) can here formally 
be derived as follows, the derivation in other physical systems (e.g., 
crack fronts) is similar. We assume the motion of the wetting line 
u: R — > R to be slow compared to the relaxation time of the water 
surface U: RxR + — > R. The system contains an energy term stemming 
from the water's surface energy, which is given (after removing the 
constant term from a completely flat surface) as 



xR+ 



•1 + \/\VU\ 2 + ldx 
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Linearizing this energy around a nearly flat state and cancelling the 
constant term we approximate this as 

I \ \VU\ 2 dx. 

</RxR+ 2 

It is well known that the infimum of this energy subject to the condition 
that U = u on the boundary of the domain is given by |[m]^i/ 2 , i.e., 
the if 1//2 -norm squared of the function setting the boundary condition. 
The variation of the if 1 / 2 -norm squared yields the term containing the 
square root of the Laplacian in equation (JTJ. The constant term F 
models the constant force with which the rough surface is pulled out of 
the water. The roughness of the surface itself acts as an obstacle to the 
evolution of the wetting line — it locally requires an additional amount 
of force to overcome a grain in the sandpaper. This is modelled by the 
heterogeneous force term /, yielding equation ([!]) as the viscous flow 
with respect to the involved force terms. An experimental example 
of this kind of system can be found in |MGR02j . They use a similar 
model, which is also proposed by [EK94, JdG84j . The non-local term 
in these articles is a mean-field version of our term. 

Another important application in which a model of the above kind 
arises is that of a crack front propagating in a rough medium. Exper- 
iments and some modeling can be found in the work of Schmittbuhl 
et. al. |SDM + 03] . The derivation of the stress intensity factor for a non- 
flat crack front (resulting in the fractional Laplacian) was first given by 
Gao and Rice [GR89J. For simulations using the model and for more 
experimental references see for example |TV04j . Nonlocal operators 
that model the interaction with elastic media also arise in models for 
dislocations |BM09l IFTMUQ] . 

In this article, we consider a specific form of the function / which is 
that of localized smooth obstacles. 

Assumption 1.1 (The random field). Fix <\>: IR 2 -»■ [0,1] € C 6 °°(M 2 ) 
with 

(i) 4>(x,y) =0 if \\(x,y)\ \ >ri, 

(ii) (f>(x,y) = 1 if ||(x,2/)|| max < r , 

for T\ > y/2r > 0. By ||-|| we denote here the Euclidean norm in R 2 , 
by I HI the maximum-norm in IR 2 . We consider f to be of the form 

J 1 1 1 1 max J 

f(x, y, u) = ^2 fk(w)<i>(x ~ x k(u),y- l/fc(w)). 

The random coefficients for the strength of obstacles fk are iid strictly 
positive random variables. The random distribution of obstacle sites 
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{(xk, Vk)}k£K is a 2- dimensional Poisson process on 1 x [?r,oo) with 
intensity A. 

Remark 1.2. Given Assumption [Li] it is clear that ([I]) admits a unique 
classical solution. 

Under these conditions we can state the main theorem of this article. 

Theorem 1.3 (Pinning of interfaces). Assume the function f is chosen 



according to Assumption 1.1 and let A := — (— A) s with s G [1/2,1). 
Then there exist a deterministic F* > and a C°° random function 
u : R x Q — y [0, oo) such that 

> (A«(w))(i) - f{x,u{x,u),u) + F* 

for almost all u £ Q and all i6l. 

Remark 1.4. By the comparison principle, of course any random field 
that can be bounded from below for a.e. w by a field of the type of 
Assumption |1.1| yields the same pinning result. 

This theorem states that there is a non-trivial pinning threshold in 
our model, since by the comparison principle any solution of the frac- 
tional diffusion problem ([I]) with F < F* and zero initial condition 
must remain below the non-negative supersolution u. Note that As- 



sumption LI ensures that an identically zero function is a stationary 
subsolution to the evolution problem ([I]) for any F > 0. Thus, for 
< F < F*, the interface becomes trapped and reaches — at least 
asymptotically — a stationary state. 

The article is organized as follows. In Section [2] we show the exis- 
tence of a non-trivial threshold for the existence of infinite percolation 
clusters which contain the graph of a function that only grows logarith- 
mically. This is a generalization of Lipschitz percolation |DDG + 10] , the 
proof for our result is inspired by |GH10j . In Section |3j using the perco- 
lation result, the supersolution is constructed. In contrast to [D DSllj . 
due to the non-local nature of the problem, a simple piecewise con- 
struction is no longer sufficient. Finally, in Section [4] we present some 
conclusions and open problems. 

2. Flat Percolation Clusters 

In this section, let |.| denote the l\— norm on IR n and denote the i th 
unit vector in IR m by (i < m). 

Theorem 2.1. Consider site percolation on Z™ +1 with n > 1: for 
given p G [0, 1], each z G Z n+1 is open with probability p and closed 
otherwise, with different sites receiving independent states. For each 
nondecreasing function H : N — > N satisfying 



PINNING IN RANDOM ELASTIC MEDIA 



5 



i) H(0) = 

ii) H(l) > 1 

iii) liminffc^oo ^ > 0, 

there exists some pu = Puij 1 ) £ (0, 1) such for each p G (pn, 1) and 
almost every realization of the site percolation model with parameter p, 
there exists a (random) function A : Z n — > N which satisfies 

a) |A(x)-A(y)| < H(\\x - y\\) for all x,y eZ" 

b) (x, A(x)) open /or every ifZ", 

Proof. It suffices to prove the theorem in case if (A; + 1) < if (/c) + 1 for 
all e No (which implies if (1) = 1). Further, we can and will assume 
without loss of generality that for all positive integers ki,...,k m , we 
have //(V /••)< V //(/.;, )• 
For j G N define 

R(j) ■= SU p{k G N : H(k) = j}. 

Now we explain what we mean by an admissible path. Let x, y G Z n+1 . 
A blocking-path from x to y is a finite sequence of distinct sites x = 
x ,Xi, ...,Xk = y in Z n+1 such that for each i = l,...,k the difference 
Xi — Xj_i takes either the value e n+ i or (z, — if(||z||)) for some z G 
Z n \{0}. A blocking-path is called an admissible path if in addition for 
each i = 1, k we have that if Xj — Xj_i = e n+ i then Xj is closed. For 
z eZ n , define 

A(z) := 1 + sup{/i G N : 
there exist x G Z n and an admissible path from (x,0) to 

If this supremum is finite for some z, then (z,A(z)) is clearly open. In 
this case, A(x) is finite for all x G Z n and A satisfies a) in the Theorem. 
All that remains to be shown is that A(0) is finite. 

By assumption, there exist C > and 7 > such that R(i) < Ce r 
for all i e N. Observe, that there exists some K such that the number 
of sites in Z n with l\— norm at most k is bounded by Kk n for all k G N. 
Fix N G N , h G N, and q := 1 - p G (0, 1). For x G Z n with ||x|| = A" 
we estimate the expected number of admissible paths from (x, 0) to 
(0, h) as follows: 

for a given such admissible path, let ki be the number of steps of the 
path containing a down-jump of size % G N. Then the expected number 
of such admissible paths which contain exactly D G N down-steps (in 
the sense that X]f=i(- r * — x i-i)~ — D) and therefore D + h up-steps is 
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at most 

D 



< a D+h e^ nD ( ( kl + "' + k ° + D + (xC n ) kl+ - +kD 
~ q h,...,k D ,D + h ) K ' 

m=0 v ' 

= 9 -.-i:((^r( m ^: k )E( fcl ,. m ,J 



m=0 
D 



< q D + h e ^ D (^KC n ) m T l+D+h 2 D ( m + D l M 



m=0 

D 



< q D+h e< nD 2 2D+h ((2KC n ) V l) D ^ ' "' : D 1 



? D+ V nD 2 2D+ ' l ((2 J ft:C n ) V 1) D 



m 

m=0 

'2D y 
D 

< q D+h ^ nD 2 2D+h ({2KC n ) V 1) D 2 2D , 

where the first two sums are extended over all ki, ko G N satisfying 
J^^Lx i/cj = £) and the fourth and sixth sums extend over all k±, ko £ 
N which in addition satisfy ^2 i=1 ki = m. Let 

(3 := 16e 7n ((2KC n ) V 1). 

Summing over D from H(N) to oo, we see that, for q/3 < 1, the 
expected number of admissible paths from (x, 0) to (0, h) is at most 
{2q) h (qf3) H{ - N \l — q(5)~ l . The total expected number of admissible 
paths starting from any point (x,0), x G 7L n and ending at (0,h) is 
bounded by 



(2q) h (l - qP)- 1 ({{KN-- 1 ) V l)(q(3) H(N 



\ U » ^ ) ■ <)\M i 

N=0 

Here, K is a constant chosen such that the number of x E Z n such 
that ||o:|| = N is bounded by KN^ 1 for all N eN. The sum is clearly 
finite provided q > is sufficiently small. Now, the first Borel-Cantelli 
Lemma implies that the largest h for which there exists an admissible 
path from some (x, 0) and ending at (0, h) is finite almost surely and 
therefore the assertion is proved. □ 
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Remark 2.2. The theorem is sharp in the sense that it becomes wrong if 
in iii) "inf " is dropped and ">" is replaced by "=" (this is a consequence 
of the second Borel-Cantelli Lemma) 

3. Construction of the supersolution 

The construction of the supersolution is performed in a series of 
steps. We first split up 1R 2 into boxes large enough so that boxes that 
contain an obstacle of a minimum strength percolate in the sense of 
Section |2j All obstacles not necessary for the percolation cluster are 
then disregarded. In each column of boxes we now have one obstacle at 
position (xi,y(xi)). Starting from a periodic supersolution (assuming 
obstacles at y — and at periodic distance in x with period larger 
than the box size), we construct a supersolution for obstacles centered 
at (xi, 0) by cutting out one period and using this function locally 
around obstacle sites. Finally, we can add a smooth function with less- 
than-linear growth (given by the percolation cluster) in order to obtain 
a supersolution that passes through the original obstacle sites. 

In this section, we make frequent use of the equivalence of the in- 
tegral representation and the Fourier representation of the fractional 
Laplacian. Furthermore, we use the symmetry of the fractional Laplace 
operator and the weak form of it by switching between applying it to a 
test function and the function itself. Further information can be found 
in |DNPVlT] . The extension problem related to fractional Laplacians 
has been treated in [US07J 

Definition 3.1. Consider thus first a,b,5,F 2 > with a > 46 and 
5 < min{l, b}. Let p := b + 5/2. Let F 1 = ^F 2 . We define 




for x E \—b — 5/2, b + 5/2] = [-p, p] 
for x E [—a, a] \ [— p, p], 



periodically extended to the real line. Now let 

g ■= vs/2 * g, 

with r/s/2 a standard-mollifiei^] with | supp r/s/2 \ = 5, i.e., radius 5/2. 

Remark 3.2. Note that the following statements hold. 

(i) g is periodic with period 2a, 

(ii) g(x) = F 2 for x E (—6, 6), 

(iii) g(x) = -F 1 for x e (-2a + b + 5,2a - b - 5) \[-b - 6,b + 5] 

(iv) g and g have vanishing averages. 



It is important that this standard mollifier used here is symmetric with respect 
to the origin. 
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Definition 3.3. In order to construct a periodic supersolution we let 
v be the modulo a constant unique and continuous periodic solution of 

Av = g 

and set v = rjs/2 * v. The constant is set so that the average of v 
vanishes. 



Remark 3.4. Note that Av = g. 

An explicit uniformly converging Fourier series representation of v 
can now easily be obtained and yields L°°-bounds on v and v as well 
as some symmetry properties. 

Lemma 3.5. We have 



i v (x 



-2« 2s ^ sin (**(& + 5/2)) ^ftg 1 



(ii) v(—x) = v(x) for all x G K, v(—x) = v(x) for all xel, o and w 
are periodic with period 2a ana 1 1> zs continuous. 

(iii) IMIoc < Plloo < ^ k2 C(2s)a 2s - 1 (b + 5/2) for s > 1/2 

(iv) Moc < (ML 

< I (Fx + F 2 ) (6 + 5/2) (2 + log (a) - log (vr (6 + 5/2))) /or s = 1/2 

ifere, £ denotes the Riemann zeta function. 

Proof. The series representation follows from a straight-forward calcu- 
lation, properties (ii)-(iii) follow directly from the representation by 
bounding the absolute value of the sine and cosine functions by unity. 
Property (iv) can be seen as follows. For s = 1/2 we have 



\v(x) 



-2a £ 



F 1 + F 2 



cos 



k^x) 

a ) 



k=l 

2 (Fx + F 2 ) 
^ 



7T^ 



fc=l 



sin (*I (6 + 512)) k2 
sin (A;* (6 + 5/2)) cos (fc^x) 



With |sin (fcf (6 + 5/2)) | < k\ {b + 5/2) and |cos(Jfe*a;)| < 1 and 
| sin (k^ (b + 5/2)) | < 1 by splitting up the sum in two parts and using 
an integral estimate for each part one gets 
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\v(x )\< { *T 2) a - 6 + 5/2 + / - (6 + 5/2) | dA:+ 



t(M-«/2) 

2 (Fx + F 2 ) 



7T (& + 5/2) U+log — - — — + 



7T 2 ^ v ' ' \ \n(b + 5/2) 

2 (Fx + F 2 



77(b+8/2) 

(6 + 5/2) (2 + log (a) - log (tt (6 + 5/2))). 

IT 

□ 

It is now necessary to establish some monotonicity properties of the 
function v. 

Lemma 3.6. The function v strictly increases on [0, a] (and thus by 
symmetry strictly decreases on [—a, 0}). 

We prove this lemma as a consequence of the following two Proposi- 
tions by showing positivity (respectively, negativity) of v' . This prop- 



erty remains valid under mollification, as shown in Lemma [3. 10 

We denote by F+ := U^l-oo ( 2ka ~ b ~ 5 A 2ka + b + 6/2) the set 
where the second derivative v" will be shown to be strictly positive and 
E~ := M. \ E + the set where the second derivative will be shown to be 
strictly negative. We also first have to show smoothness of v on the 
union of those two sets. 

Proposition 3.7. Let p 6 (0,1). We have that (—A) p g(x) (given by 
its integral representation) exists for all x G E + U E~ and one has 



■A)^(x) 



> o x e F^ 
< o x e E~ 



Proof. Take x e E + , then one has (with some C > 0) 



:-A) P 3(x)=CP.V. / ? + %' dy 

\x — y\ 



g(x)-g (y) 

\x — y\ 

^ RV - / | Fl + l f+2 P XE + (y) dy>0 
Jr \x — y\ 
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where xe+ is the characteristic function of E + . For x G E~ the same 
calculation gives 



-Ay g (x)=cp.v. / TT+^dy 

7r \x — y\ 



x-y\ 
x-y\ 



-CP-y- I , Fl + n+ 2p XE- (y) dy<0 



□ 



Proposition 3.8. For all x G E + U E the function v is twice differ- 
entiable with 

v" (x) = Ad (x) = (-A) 1 " 5 (- (-A) s d (x)) = - (-A) 1 " 5 g (x) . 

Proof. Let B := (-A) 1_s 

For B given by the integral representation Bg is continous on (— p, p) 
and [—a, a] \ [—p, p\. For [— f3] C (— p, p) take e > so small that 
[-(3-2e,(3 + 2e]c(-p,p). 

Take a standard mollifier i] t with radius e. Define v e := v * i] t and 
take a test function tp G with support suppip C [— /3, /?]. Then one 
has for x G [— /3, /?] 

v'^x) = Av € (x) = -B(Av t (x)) — - B(A(v * r) e (x))) 
= -B((Av) *ri e (x)) = -B(g*r] e (x)) 

As g * i] t is smooth and bounded one has 

(-B(g*r) e ), ip) = (g*r) e , -Bip) = (g, — (Btp) * i] t ) 

= (g, - B(^j*r] e )} 

Using this together with the integral representation of B one gets up 
to some constant 

2{g,B(i(;*r ]e )) 

, . ij*r] t (x) - ip*ri e (y) \ 

(J( } ' ' — | X _ y |i+2(M — d yj dx 

g(y) I / — . ^2^) — dx d y 

r \x — y\ v ' J 

t^z^KW * Ve{x) -y>* Ve{y))dxdy 

\x - y\ 

-^^ip * r] e {x)dydx - / / -^^ip * r] t {y)dxdy 

\x - y\ Jr Jk \x - y\ 
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( f g(x)-g(y) A , 

-/3-e \J]R \X — y\ 



g{x) - g 



i>*r)e(y)dy 



\ I I |3-2s / 

-p-e \Jr \x — y\ J 
= 2(Bg,^* Ve ) 

From the choice of [—/3, (3], e as well as suppip C [— /3, (3} one has 
(Bg, ip *f] e ) = / Bg(x) / ip(y)i] t (x - y)dydx 

J-P~£ JP 

rP+e 

B g(x)ip(y)r) e (x - y)dydx 



P-eJ-P 
P+e rP 



j B g{x)ip(y)rj e (y - x)dydx 

■P-e J-P 
rP rP+t 

I 4>(y) Bg(x)r] e (y - x)dxdy = (Bg*r) e , ip) 

J-P J~B-e 



-P 

Finally this shows that 

K( x ) = - B (9 * Ve(x)) = (-Bg)* rj e (x) 

for all x G [-/3, 0]. 

In the same way one can show that 

K( x ) = ~ B (9 * Ve(x)) = (- B g) * ^(x) 

for all x G [-£, /3] \ [-p, p] with [-£, f3] \ [-p, p\ C [-a, a] \ [-p, p] 
From the continuity of — B g on [— /3, /3] one gets that v" = (— B g)*rj e 

converges uniformly to — Bg on [— /3, (3]. 
Furthermore one has 

(i) v e — > 5 uniformly on [— /3, /3] 

(ii) v G if^H?, /?]) and therefore v" = Av G iJ" 1 ([-£, /?]) 

Then one has for any test function ip G with support in [—/3, 0\ 

(v, i;") = [im (v e , r) = [im = (~ B & V) 

This shows that {;" = — Bg in /3, /?]) and in the sense of 

distributions. As — B g is continuous on [—(3, (3} this proves that v G 
C 2 ([-(3, 0\). 

As [— /3, (3} C (— p, p) was arbitrary one gets 



v"(x) = -B~g(x) = -(-A) 1 - s ~g(x) 
for all x G (— p, p). 
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In the same way one shows the assertion for x G [—a, a] \ [— p, p] and 
from periodicity the statement follows. □ 

Remark 3.9. For symmetry reasons we have v'(0) = v'(a) = 0. 



Proof of Lemma 3.6 , According to Proposition v" exists on E + U 
E~ and is strictly positive on E + and strictly negative on E~ . Be- 
cause of v' (0) = = v' (a) we get for x G [0, a] \ {b + 5/2} by the 
fundamental theorem of calculus 



v' (x) 



J^(-A) 1 - S g(y)dy xE [0,6 + 5/2) 
-f«(-A) l - s ~g(y)dy xe(b + 5/2,a] 



> 



where in the last step Proposition |3.7| was used. □ 

Lemma 3.10. The function v strictly increases on [0, a] (and therefore 
by symmetry strictly decreases on [—a, 0}). 

Proof. Note first that v' G Lf oc as a weak derivative. Since 

V = —{y * r) S ) = — V * 1]s 

ax \ ax J 

one gets 

v ' i x ) = ( * % ) ^ > 



for all x G (0, a) according to Lemma 3.6 and because the convolution 
preserves positivity in this case due to the symmetry properties of v 
and the mollifier. 

□ 

We now split 1R 2 into boxes large enough so that the percolation the- 
orem from Section [2] can be applied to boxes that contain an obstacle. 

Definition 3.11. For k G Z, j G N and I, d, h > 0, I > 2r x let 



Qk 




-d) - 


■ l/2,k(l + d) + 1/2} 


Qk 


:=[k(U 


-d) - 


-l/2 + r 1 ,k(l + d)+l/2-r l ] 


Qkj 


■=Qk x 


[(;- 


l)h,jh] 



The following is a direct result from Section [2] and Assumption 1.1 



Proposition 3.12. Let < a < 1, H(k) := [k a \ (i.e., the integer 
floor of k a ) and let V := (7 — 2r\)h > 0, q > such that 

l-exp{-\V-P{f 1 >q}}>p H (l) 
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from Theorem 2J_, Then, almost surely, there exist a random function 
A: Z ->■ N with \A(x) - A(y)\ < H(\x-y\) for all x,y E Z and a 
mapping I:7L^-K with 

(%i(k),yi(k)) e Qk,A(k), fi(k) > <? 

/or a// fc£Z. In £/ie following we denote by I the set /(Z). 
Definition 3.13. Let now d > I, 2a > d + 21 and define for isl 

/ \ f u(a; — Xj) for a; G [sCj — Z — d/2, X{ + Z + d/2] 
^ ' ^ +oo otherwise. 

:= f (x — Xj) 

Wflat(^) := mva.Vi(x), 



where the 2a-periodic function v. R — )■ R is given in Definition 3.3| 



Remark 3.14. Note that due t o the monotonicity and periodicity prop- 
erties^] of v from Lemma 3.10 we have that on any interval of the form 



(6, b + 2a] two functions Vi, Vj, i,jEl intersect each other exactly twice 
or they are identical. The points of intersection have distance a. 

Proposition 3.15. With the definitions above we have that 

(i) M flat is bounded and continuous, 

(ii) given Xi e Q k , we have that UQ a t(x) = Vi(x) for all x E Q k - 

Proof, (i) is obvious by construction. In order to see (ii), note that the 
spacing d between two boxes Qj is larger than the length Z of a box 
and the period of the function v is larger than d/2 + Z. □ 

In the following, we prove that the function Wfl at constructed above is 
a supersolution to a modified problem where the obstacles are extended 
from — oo to +oo in the y-direction. We thus fix £ E R and calculate 
the effect of the fractional Laplacian A on u^ at evaluated at £. In the 
case that £ is a point where Wfl at is smooth, i.e., not a point where 



the minimizing Vi in Definition 3.13 changes, we can directly apply the 



integral representation of A. The points of discontinuity (of the first 
derivative) of Wfl at will have to be smoothed in order to construct a C°° 
supersolution. 

Definition 3.16. Given £ <E R, let i E I such that Ufl a t(£) = v i (0 = 



u 



(£). If Vj(£) = Vi(£) = Mfl at (0 we take i = max{i,j}. Further- 
more, we recursively define the points of intersection of the periodic 



2 The function v admits one maximum and one minimum in each period and it 
is strictly monotone in between. 
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supersolution w io with u& a t- Let 

• / ^ n =i ^ n u io> M flat on (£ - y - k, £ - y) 

ai := mm < y > : =k > with " ;1 , C 

[ ~ u io < M flat on (£ - y, £ - y + K) 

u • / ^ n ^ n < M flat on (£ - y - k, £ - y) 

6i :=mm < y > ai : 3k > with Yl ^ , ( 

[ ~~ Mi > M flat on (£ - y, £ - y + k) 

a k+1 : = min (y > b k : 3k > with > Mflat ° n % ~ V ~ «' * ~ V \ 
+ \ y ~ Ui < Mflat on (f - y, «e - y + k) 

h — min/i, >n • 3k >> with < M flat On (f - y - K, £ - J/) 

.-mm |y > a fc+1 . 3k > with ^ > ^ on (£ y, £ y + k) 

Define a*., b k in the same way by substituting — y by +y 
Lemma 3.17. We have 

b k ~ a k > a owe? Ofc+i — b k < a 

for all fceN. 

Proof. Take G N arbitarily and ii,*2 G ^ such that — a k ) = 
i*flat(£ — Qfc), v i2 (£ — = Mflat(£ — (it is clear by constrcution that 
i\ and i 2 are unique). 

Because of Remark |3. 14 and the construction of u^ t one has 

Wflat < «ii on [C - Qfc - °> C - fl fc] • 
The definition of a k and Remark |3. 14 yield 

Ui (£- a k) = «ii and «i ( a; ) > ^u^) for ^ £ (£-afc-a,£-a fe )- 

Altogether one has 

M io (x) > Mi^x) > M flat (x) for all x G (£ - a k - a, £ - a k ). 
Furthermore 

Mflat < Ui °n [6 - &fc, 6 - °fc] 
and by the choice of and 6fc there exists no larger interval J D 
[£ ~ b k ,!; — a k ] with «fl at < u io on J. Using this one gets 

(6 - a k - a, £ - ajfe) C [£ - 6 fc , £ - a fe] 

and therefore 

b k - a k > a. 

The other inequality is shown by an explicit calculation. By construc- 
tion of M flat it is obvious that for I 3 i 3 :— min{i G I\i < i 2 } one 
has v i3 (£, - a k+1 ) = M flat (^ - a k+1 ). Define r := 2 (x i2 - - 6*)) and 
z := — r - Then by the periodicity property of u io and u i2 it is clear 
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Figure 2. Illustration to estimate the distances be- 
tween intersection points of local supersolutions with the 
periodic supersolution Vi centered around Xq. 

that Zp is a minimum of Ui . Furthermore from the proof of Propo- 



sition 



3.15 



one knows that Vi 3 and fj 2 intersect in ' 3 ' 2 . The same 



argument also shows for the intersection of Ui and v l3 

X;.. -\- Zn Xt- 4- Xa. — V 

£ — O-k+1 



2 2 2 2 

By the choice of the Qu it follows \xi 2 — Xi 3 \ = Xi 2 — Xi 3 < d + 21 — 2t\ 

and therefore — ■ < d/2 + I — r\ < a. Putting everything together 
one gets 

flfc+i - Ofc = (£ - 0*) - (? - = g 2 = 2 ° 

See Figure [2] for an illustration. 

□ 

Lemma 3.18. Assume that £ is not a point of discontinuity of the first 
derivative o/«fl at . VKe then have 

\y -CI •'K 12/ — CI 

Proof. We only consider the part of the integral to the left of £. With 
the definition of a^fefc, the estimate for the other part of the integral 
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follows accordingly. Using the same notation as in Lemma 3.17 



one 



has itflat < Ui on [£ — £ — a*] and so the definition of Wfl at together 
with Lemma 3.17 yields 



ttflat (y) - u io (y)dy 



u i2 (y) -u iQ {y)dy, 



u i2 (y) -u io (y)dy < 0. 



-o-k 



Un{y) -u io (y)dy 



< min 

< / 

In the same way using u io < u^t on [£ — Ofc+i, £ — bp.] one gets 



< 



Mflat (y) - U io (y)dy 



< min 



£-a fc+1 +a 



(y) -u io (y)dy, 



< 



Ui 2 



t-b k 



u i2 (y) -u io (y)dy 



f-6 fc -a 



Using this one can split up J^_^ +i ""j^i 1 "' 2 "^ d y * n ^ wo P ar ts where 
the integrand is negative and positive, respectively. Noticing that 

^ ^ on K _ a fc+i' £ _ and < | y _gji+^ on [£ - 6*, £ - a*] 

one gets 

/"^ ttflat(y) -Wio(^) 



I?/ — ei 



l+2s 



dy < 



+ 



< 



bl +2s 



bl +2s 



Z-b k 

Mflat (y) - u io (y)dy 
Mflat (y) - u io (y)dy 
w i2 (y) -u io (y)dy = 0, 



where in the last step the peridicity of Wj and w; 2 was used. See Figure [3 
for an illustration. Inserting a zero in the form — m (£) + Wi (£) one 

gets 



> 



r+2^ d ?/ 



£— a k +i 



|y-CI J 
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Figure 3. 



Using now Ui (£) = Ufl at (£) and summing up for all k in N it follows 
that 

r( - ai u Rat (y)- Ufiat (0 A u i0 (y)- Uio (Q 

Furthermore an analogous calculation shows that 



Wflat (V) - "flat (0 



'n+ai \y-£\ 1+2s " Jn+ai Is/ -61 

which together with Wfl at = Ui on (£ — ax, £ + 5i) yields 



1+2^ — d 2/' 



"flat (j/) - "flat (0 ^ < ^ ^0 iv) ~ ^0 (0 dy 



is/ — ei 



l+2s 



□ 



Next we show that the function Ufl at is a weak supersolution to the 
above mentioned modified problem. 

Lemma 3.19. Let < e < ^, 6 + 5 < r Q - 2e, F 2 < q and let 

<7fi a t : M — > [0, oo) &e gzuen as 

, v f 9 /or x e Uie/fo ~ r o + l e ' ^ + r o - |e] 
^ Hat ^ I for x <£ {J ieI (xi - r + §e, x< + r - |e) . 
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Then for all ip E C C °°(M), ip > we have 

(ttflat, Aip) + (-g m + F,ip) : = / M flat (Aip) dx+ / (-o flat +F) ip dx < 0, 

</r </r 

/or any F < min{g — F 2 , Fx}. 

Proof. If the support of ip does not contain any points of discontinuity 



of 

u fiat' the statement is clear from Lemma |3.18| by noting that the 
integral representation of Au^t is well defined and finite on the whole 
support of ij} and one can thus apply the fractional Laplacian directly 
to Mflat- A calculation then yields 

^flat (y) - Mflat (x) 



\y - x 



l+2s 



dy - Oflat (x) 



< 



^ / , n+2i dy- o flat (x) 

Jr \y — x\ 

F 2 - Oflat (x) for x E [x i{x) - b-5, x i(x) +b + S] 
-Fx - Ofl a t (x) for x E [-a, a] \ [x i{x) -b-5, x^ x) + b + 5] 

< -min{g - F 2 , Fi} . 

Using a partition of unity, one can isolate any points of discontinuity. It 
is possible to locally split Wfl at into a piecewise affine function with the 
same jump and a C 1 fl Hf oc -i unction with bounded second derivatives. 
Noting that the first derivative of Ufl a t always admits a negative jump, 
the integral operator A applied to those two parts yields a strongly 
negative term near the jump for the piecewise affine function and a 
bounded term for the C 1 fl Hf oc -f unction. The statement in the lemma 
is thus proved for all test functions ip. □ 

It is now possible to mollify the function «fl a t by a standard mollifier 
of radius e to obtain a smooth classical supersolution to the modified 
(flat) problem. 

Corollary 3.20. Let M smooth := Ve * "flat e C°° and g smoot h ■= ij e * Stoat- 
We then have 

A ^smooth — ^smooth + F < 

for all F < min{g — F 2 , Fi}. 

Proposition 3.21. Let h > 0, d > 0, I > 0, s E [1/2, 1). For A : Z -> 

K such that |A (zi) — A (z 2 )| < 2h\z\ — z 2 | a with < a < 1 there exist 
a smooth function w s tcp : M — > M and constants Cq, C\ and C 2; that 
only depend on s and a such that 

(i) M step (a;) = A(k) for any x E Q k , k E Z 
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ii) HdXtcpIL < CoJ> 

fl ^ ( ~" 2 (d/2+l/2f 



(iii) | (- A)' u stcp (x) | < d W 2 +V 2 ) 2 ' 2s h + C 2 



\2s 



Proof. Parts (i)-(ii) are immediate by mollifying a piecewise constant 
function. Part (iii) can be seen as follows. Without loss of generality 
one can assume x G [—d/2 — 1/2, d/2 + 1/2]. 

Let II := (-3 (d/2 + 1/2), 3 (d/2 + 1/2)), then one has from the as- 
sumptions on A and from (p} for all y G M\H 



(2) l^stcp (x) - u stcp (y)\ < 6h 



f - y\ 

(l + d) c 
because 



i < > 



\u 



x-y\ 



step 



(x) - u stcp (y)\ < 2h a + Ah 



but as y G one has > 1 and therefore 4/t < . 

Using that u step is smooth and grows less than linear, (— A) s it step 
can be a represented as 



a \s / v W «sto P (a; - y) + u stcp (x + y)- 2u stcp (x) 
-A) u step (x) = - / dy 

1 \y\ 



by applying some standard estimates [D NPVllj one gets 



|(-A) s n stcp (s)| = - 



u stcp (x - y) + u step (x + y ) - 2u stcp (x) 



I |l+2s 

I y I 



r n \y\ 

1 f |^ s tep (g - y) ~ U stcp (x) | + |-» step (x + y) - U siep (x) \ 

oi i ii+2s y~ 

/./ro\n \y\ 



20 PATRICK W. DONDL, MICHAEL SCHEUTZOW, AND SEBASTIAN THROM 

Putting in estimate Q for M step and calculating the resulting integals 
it follows that 

|(-A) s ?i ste p(:r)| 

i-3(d/2+J/2) 1 6h f 1 

< ll^^stcpIL o2-2s / . /2 , , /o^- 2s _|_ 12/1 /" X 7/- 1 - 2s+a d77 

< ll^stcpIL „ a _a. , , //9 ,2-2s , 12/1 o2s - a (rf/2 + V2)- 2S+a 

" 2-2 S d W2 + 1/2) + 2s _^ (d + i) « 

* ^ 2 + + a^T?* ^ 2 + ^ 2s 

where in the last step ll^Wsteplloo < Co J? was used. For 

q2— 2s g2s— « 



Ci = Cn and C 2 = 

2 -2s 2s-a 2° 

one obtains the estimate. □ 

Lemma 3.22. Let s > 1/2 and take := -liC (2s) and < C s < 1, 
C a >5. 



Take q > and > as in Proposition 3.12, Choose < _F 2 < q and 



take F\ > as in Definition 3.1 Choose now I > such that 

I > max I Ar u V(2S) , (12 + C7 a ) Vr f^ 

1 + 2F 2 r n + 12F 2 (d + C 2 ) V^pCg 1 

' (M!r^) 1/(2S_1) 

. Take d = l, ll = I + l < a < CJ, b = § (CooF f a ° 2s+ro) , < 5 < C 5 6 
and ft, = ^z^r- Choose < e < — swc/i inai 6 + 0" < r — 2e (possible 
due to item M)j. Taen we /iai>e 



(i) p = 6 + 5/2<f <£, 

(ii) (Ci + C 2 )^7zk<^ 



3 ^18' 

L_ < 

Z 2s l-2n ^ 2 

J L_ <• J_ 

Z 2 " l-2n - 12 (C*ooC 2s F 2 Z 2s +ro) 



(hi) ( Cl + C 2 )V 1 hj^- i <^ a 

(iv) ||n|| oo <^±^C(2 S )a 2 - 1 (o + 5/2) 
= C oc (F 1 + F 2 )a 2s - 1 (o + 5/2) < 



6 (CocC2"F 2 «2=+r ) ' 
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Proof. (i) The choice of b, 5, a and / yields 

h j- A /o <r oh 1 ar ° / 1 r ° 
b + 5/2 < 2b = - < 



< 



3 CooFaa 2 ' + r 3 C^a 2 '- 1 
1 r 



< 



r 1 



3 C^F, (I) 25 " 1 

To 

3 



from the conditions on r , ri and the choice of / one also 

r T\ I 2 a a 
~3 < 3" < 12 ~ 3 12 ~~ 18' 

(ii) From the conditions on Z we get 

'2(d + C 2 )y 



Z > 

z 2s > 



2n (q - F 2 



2s ^ 2(C 1 + C 2 )V 



2n (g - F 2 ) 
Z 2s (Z - 2n) > 2r x Z 2s > 



g-F 2 

where in the last step Z > Ar 1 was used. This gives 

q ~ F2 >(c 1 + c 2 )v- 1 1 



2 v P s l-2n 
(iii) From the condition 

z 1 + 2F 2 r n + 12F 2 (Ci + C 2 ) VC^C%_ 
F 2 r 

on Z we get 

(3) r F 2 l - 2F 2 r n - 12F 2 (d + C 2 ) VCJ^ > 1. 

By rearranging some terms we get from the condition 

(12(Ci + C 2 )Wo)£ <Z 
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on I that 

12 (d + C 2 ) Vr < l 2s 

< l 2s (r F 2 l - 2F 2 r ori - 12F 2 (d + C 2 ) VC^C 2 /) 
11 

12 (d + d) W < r F 2 / 1+2s - 2F 2 (r r a + 6 (d + d) VC^Cf ) / 2s 
/ 1+2s r F 2 - 2r nF 2 l 2s > 12 (d + d) VO^Cf Z 2i + 12 (Ci + Vr 



Z 2 * Z - 2n 12 (CooC^P* + r ) 



fiv) From Definition 3.1 we know 



i*i + F 2 = I — P — + 1)F 2 = -^—F 2 . 



a — p J a — p 

Furthermore from Lemma 13.51 one has 

NL<^i^c(2 3 ) a 2 - 1 (i, + «) 

Putting this together and using item dij) one gets 



< Coo^^a 28 - 1 (6 + 5/2) 

18° 

< ^dc^a 2 *^ 



17 00 G^F^+ro) 

18 „ rn ,28-ln flr 



< -dci^ 5 -^- 



17 



17 °" 6doF 2 a 2s 
6 



(v) Because F\ is as in Definition 3.1 one has 

—Fi = -^F 2 = - " + ^ 2 F 2 < - J^F 2 < -<+ 
a — p a— (6 + 5/2) a — 6 a 

1 r a 1 



ae^F^ + ro) 6(dcd 2 ^ 2 / 2s + r ) 

□ 
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Take now u := Wsmooth+Wstep- Choose the parameters as in Lemma 3.22 
and 

then one has u > as just seen and we can now give the 



Proof of Theorem \1.5\ for s > 1/2. Let the parameters be as in Lemma 3.22 



and take M ste p (x^ = jji for all i E I which, due to Propositions 3.12 



and 3.21, is (almost surely) possible. From the choice of ga at and r/ e we 
have 



~y]fi M 4>{?- X h ^smooth (x) + U stcp (x) - JJi) < -5-flat * T] e (x) 



iei 



Using this we have for F < F* 
Au (x) — f (x, u (x, u) , to) + F 

< A M smooth (x) - fi (u) <f){x- Xi, Msmooth (x) + U stcp (x) - Jji) 



+ F + Au stcp (x) 
< A u smooth (x) - g-flat *r] t (x) + F + A u stcp (x) 



With the results of Corollary |3.20| and Propositions |3.12| and |3.21| it 
follows that 



Au(x) — f (x, u (x, u) , to) + F 



< -mm{q-F 2 , Fx} + F + C x 



{d/2 + 1/2) 



2-2s 



(P 



-h + C 2 



{d/2 + 1/2) 



2 s 



< - min {q-F 2 ,F 1 } + F+ [C x + C 2 ) 



h 

J2s 



where i n the last step d = I was used. Applying now the estimates of 



Lemma 



3.22 



and h 



A u (x) — f (x, u (x, u) , u) + F 

r 



< — min < q — F 2 



6 (CcoC^FzP* + r ) 
r 



F 2 }+F + (C 1 + C 2 ) V 



1 



l 2s I - 2n 



= 0, 



which finally concludes the proof for s > 1/2. 



□ 
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For the case s = 1/2 some changes in the choice of parameters have 
to be performed due to the worse L°° estimate on v in Lemma 3J3 in 
this case. 



,3 



Lemma 3.23. Let C p = \J — - - t s 3 _ n and choose < Cg < 1, 



C a >5. 

Take q > and V > as in Proposition 3.12. Let < F 2 < q and 



choose F% > as in Lemma 3.5[ Choose now I > such that 



n (q-F 2 ) J 7 V F 2 C P 



2 



12 {Cl + C2)Ca V + 2 ri 
r 



Take d = I and |Z = |+Z < a < C a l, b = | min 



1 r 



< 5 < Csb und h = l X Lri ■ Finally choose < e < ^ such that 
b < r — 2e (possible due to item |wp/ T/ien we /icwe 

(i) Z > 4r 1; 

(ii) p = b + 5/2<b + 6 < f < § , 
(hi) (Cx + C 2 )^^<^k ; 

(iv) (C 1 + (7 2 )y |I ^ F 4 T < |F 2 min|i 

(v) iklL < ^, 



1 rp 



(vi) -F x < -F 2 min < \ 



r l L_ jro_ 



2 A/48e2(^) 3 C3' ^ 3/2 ' 6C a Z f 



Proof, (i) This is clear because I > ( ^g 1 ^ ) +4ri and ( 2y( g+ C2) ) > 


(ii) From the choice of 5 and 6 it is clear that 6 + 5 < 2b < ^. From 
the conditions on r , r\ and the choice of I we further have 
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(iii) From the condition on Z some calculation gives 

2 (Ci + C 2 ) 1/ 



/ > 



Z (/ - 2ri) > 2r x Z > 



2n (g - F 2 ) 

2 (Ci + C 2 ) V 



2 



q - F2 >(c 1 + c 2 )v 1 1 



2 v 1 ^ ll-2r x 
(iv) From the condition 

2 



on Z one obtains 

^— : < Z - 4n < / - 4n + — 



. ,,,, , <Z-4r x <Z-4r : + 
Z 2 -4nZ + 4r 2 (Z-2r x ; 



Z Z 
Therefore taking the square root and expanding by Z one gets 

2V (d + C 2 ) I - 2n _ Z (Z - 2n) 



F 2 C P Z 1 / 2 Z 3 /2 ' 

which finally, after rearranging and putting in C p , yields 



11 1 „ „ 1 1 „ 1 vrr 



{Cl + ° 2) V 1T^[ < " 2 F2 2 \/ 48e 2 (ff) 3 C a 3^' 

The second part simply follows by rearranging the condition 

l>12 {Cl + C2)Ca V + 2 ri 
r 



on Z such that 



and finally 



l-2r 1 >l2 {Cl + C2)Ca V 
r 



(v) From Lemma 3.5 we have 



ML < 2^^p (2 + log (a) - log (vrp)) . 
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The choice of F\ from Definition 3.1 together with item ^ yields 

f 1 + f 2 =(^ + i)f 2 = —f 2 <-^f 2 = ^-f 2 . 

\a — p ) a — p j^a 17 
Taking both together one gets 

Halloo < hn^r )p ( 2 + lo s (°) - log M) • 



* 17tt r 

From the choice of b and 5 we further get 



p < 2b < min 
and therefore by squaring 

p 2 < min 
so in particular 



vrrjj 1 r 



48e 2 (^) 3 v^' 3 



3 „ ' 



48e 2 (S) a 9 



3 



2 < 7Fr 



Using this by rearranging and adding additional terms we get 



ae 2 1 Tq 



2 

?"o , 1 / r \ , 1 / r 



2(^)p 2 V2(^) P ; 6 \j! 

where in the last step the standard estimate for the exponential 
function was applied. This shows that we have 

2 / ( r 



2 + log (a) < log (vrp) + 



Applying the logarithm gives 

r 



2(ft)P 
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and this finally yields 
'36F 2 



17 71 



p (2 + log (a) -log(vrp)) < 



r-2 



(vi) The choice of F\ from Definiton 3.1 and the definition of p give 

b 



-F x = - 



P 



a — p 

F 2 . 
— mm 



-Fo < 



r IF 2 . 

-r 2 = mm 

a 2 a 



Tin 




TfT'r 



1 r 



48e 2 (^) 3 ^' 3 



48e 2 (36F±\ 3 a 3 / 2 ' 3a 
As a < C,,l we have — - < — 7^-7 and therefore 

— a a — qj 

IF 2 . 



-F x < 




r 



48e 2 (^) 3 ^' 3 
1 



r 



48e 2 (ff ) 3 C7 3 ^ 3/2 ' 3CJ 



□ 



Proof of Theorem 1.3 for s = 1/2. Choose the parameters as in Lemma 3.23 
and define u := u sm ooth + «step- Then Lemma 3.23 shows u > and for 

1 



F* = - min < g — F 2 , F 2 min 



1 



TTrn 



1 



r 



2y 48e 2(^) 3 c <3Z 3 / 2 ' 6CJ 

using the estimates of Corallary 3.20 , Proposition |3 . 2 1 1 and Lemma 3.23 
we can show analogously to the proof of the case s > 1/2 that 

A u (x) — f (x, u (x, u) , u) + F < 

for all F < F*. □ 



4. Conclusions 

In this article we have shown existence of a non-trivial pinning thresh- 
old for interfaces in elastic media with local obstacles. Models of the 
kind discussed frequently arise in physics, for example in the propaga- 
tion of crack fronts in heterogeneous media. Assuming free propagation 
of such an interface for large enough driving force (which is trivial to 
obtain under some conditions on the heterogeneity), we have shown the 
transition of a microscopically viscous kinetic relation (force=velocity) 
for interfaces in elastic media with random obstacles to a stick-slip be- 
havior on larger scales. The construction of the supersolution has been 
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constrained to the 1+1 dimensional case, i.e., that of a 1-dimensional 
interface propagating in a 2-dimensional plane. In many cases, this is 
the physically relevant situation. The n-dimensional case is still open 
due to technical difficulties concerning mostly the compensation of er- 
rors arising from modifying a periodic solution. 

Furthermore, we have shown a percolation result, namely a non- 
trivial percolation threshold for the existence of an infinite cluster in 
next-nearest neighbor site percolation that is the graph of an only log- 
arithmically growing function. 
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